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We report a new locality for Pedicularis spicata Pall. (Orobanchaceae) in the Hakusan Mountains, Hon¬ 
shu, Japan. Although the populations of Hakusan Mountains have been described as P. verticillata , we 
showed that they are P. spicata. Pedicularis spicata is identical with the Hakusan populations in trnK 
sequences of chloroplast DNA (cpDNA), whereas internal transcribed spacer (ITS) phylogeny unites a 
clade of P. spicata and the Hakusan populations at 94% bootstrap support. Morphological analyses fur¬ 
ther suggest that the population at Hakusan Mountains should be included in P. spicata rather than P. 
verticillata. Pedicularis spicata occurs mainly in northeastern Asia. This study shows a disjunct pattern 
for P. spicata between the Asian mainland and Honshu, Japan. 
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The flora of the Japanese archipelago has been 
relatively well investigated with the results sum¬ 
marized in the floristic literature (Ohwi 1965, 
Iwatsuki et al. 1993, 1995a, 1995b, 1999, 2001, 
2006). Many new taxa and geographic occur¬ 
rence records, however, continue to be reported 
(Setoguchi & Maeda 2010, Oda & Nagamasu 
2011, Kadota 2012, Suetsugu et al. 2012, Suga- 
wara 2012, Fujii et al. 2013). Such findings are 
important, because they can prompt the reexami¬ 
nation of previous taxonomic assignments and 
delimitation of taxonomic boundaries. Further¬ 


more, improvements to the quality and accuracy 
of geographic distributions are important for un¬ 
derstanding the biogeographic origins of the Jap¬ 
anese flora (Murata 2000). 

Pedicularis verticillata L. (Orobanchaceae) is 
an annual alpine herb that has been classified into 
ser. Verticillatae Maxim, of sect. Sigmantha H. L. 
Li (Yamazaki 1993). This species is widely dis¬ 
tributed across the high latitudes of the northern 
hemisphere and thus is a circumpolar arctic-al¬ 
pine plant. In Japan, it shows a restricted distribu¬ 
tion pattern and has been recorded in the Taisetsu 
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Mountains (Mtns.) in Hokkaido and on several 
high mountains in central Honshu (Mt. Shirou- 
ma, Mt. Shibutsu, Yatsugadake Mtns., Akaishi 
Mtns., and others) (Yamazaki 1993, Shimizu 
1997). The occurrence of Pedicularis verticillata 
on the Hakusan Mtns., Ishikawa Prefecture, cen¬ 
tral Honshu, Japan, has also been reported in the 
literature (Hashimoto & Satomi 1969, Satomi 
1983). In a recent detailed inventory of the moun¬ 
tains, P. verticillata was reported along with the 
description of a voucher specimen (Hakusan Na¬ 
ture Conservation Center 1995). The populations 
occur in the subalpine-alpine zone at an eleva¬ 
tion of 2,000 m or more. It was observed, how¬ 
ever, that the morphological characteristics of the 
plants (size and inflorescence traits) slightly dif¬ 
fered from P. verticillata on other mountains in 
Japan (N. Shirai & T. Shimizu unpublished). Fujii 
et al. (1996) examined variation in chloroplast 
DNA (cpDNA) among the populations of P. verti¬ 
cillata in Japan, including those on Hakusan 
Mtns., and observed a large sequence differentia¬ 
tion between the Hakusan accessions and those 
from other mountains. Based on the results, it 
was inferred that the Hakusan populations might 
not be P. verticillata. 

To clarify the identity of the Hakusan popula¬ 
tions, we reexamined them using the taxonomic 
literature and molecular phylogenetic analysis. 
The examination suggested that the Hakusan 
populations were Pedicularis spicata Pall (Fig. 
1), a species, like P. verticillata, classified in ser. 
Verticillatae of sect. Sigmantha (Yamazaki 1993). 
Pedicularis spicata is confined to eastern Asia, 
ranging from Hokkaido to northern Korea, north¬ 
eastern China, northern China and Mongolia to 
eastern Siberia (Ivanina 1991, Yamazaki 1993, 
Yang et al. 1998, Grubov 2001) (Fig. 2). On the 
Japanese Islands, P. spicata has been reported 
only in the Tokachi district of Hokkaido. It was 
suggested that it might have been introduced 
from northeastern China (Yokoyama 1951, Takita 
2001). Until now, P. spicata has not been reported 
on the Honshu Island, Japan. 

To determine the identity of the Hakusan 
plants, we performed molecular phylogenetic 
analyses and morphological comparisons with 


closely related species of Pedicularis ser. Verti¬ 
cillatae ( P. spicata, P. verticillata, and/', refracta 
(Maxim.) Maxim.) in Japan. 

Materials and Methods 

Materials 

For molecular phylogenetic analyses, DNA 
sequences of the nuclear internal transcribed 
spacer (ITS) region and the plastid trnK 5’-termi- 
nal region (including a part of matK ) were deter¬ 
mined for three accessions from the Hakusan 
populations, one accession of Pedicularis spica¬ 
ta, five accessions of P. verticillata, and three ac¬ 
cessions of P. refracta, excluding the accessions 
from Mt. Kita and Mt. Tolbachik (Table 1). All 
sequence data were deposited in the DDBJ/ 
EMBL/GenBank databases under accession 
numbers AB754816-AB754828. The new se¬ 
quence data, including the ITS and trnK sequenc¬ 
es, were assembled into alignments of eight spe¬ 
cies (15 accessions) previously published by Ree 
(2005) and Yang & Wang (2007). Of these spe¬ 
cies, P. lineata Franch. ex Maxim, and P. szetsch- 
uanica Maxim, form a strongly supported clade 
with P. verticillata, whereas P. anas Maxim., P. 
megalantha D. Don and P. longiflora Rudolph 
represent outgroups (Ree 2005). Voucher speci¬ 
mens newly used in the present study were depos¬ 
ited in the herbarium of the Faculty of Science, 
Kumamoto University (KUMA); Makino Her¬ 
barium of Tokyo Metropolitan University 
(MAK); and the herbarium of the Arnold Arbore¬ 
tum, Harvard University (A). 

For morphological analyses, we examined 
herbarium specimens at the University of Tokyo 
(Tl), Kanazawa University (KANA), and the Bo¬ 
tanical Garden of Hokkaido University (SAPT). 
In addition, we used dried specimens prepared 
from our own field collections of the Hakusan 
populations and of Pedicularis refracta (from 
Takamori-cho, Kumamoto Prefi, Kyushu, Japan). 
In total, we examined 20 specimens of the Haku¬ 
san populations, 127 of P. verticillata, 44 of P. 
spicata, and 23 of P. refracta. Voucher specimens 
for our collections were deposited in KUMA. 
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Fig. 1. Pedicularis spicata in the Hakusan Mtns., Honshu, Japan. A. Habit (Mt. Shaka-Mt. Nanakura, August 
14, 2006). B. Close-up of flower (Mt. Shaka-Mt. Nanakura, August 5, 2010). Scale: 5 mm. 



Fig. 2. Distribution map of Pedicularis spicata. The distribution range on the Asian continent roughly is shown 
on the basis of Grubov (2001), Ivanina (1991), Yamazaki (1993), and Yang et al. (1998). In Japan, this spe¬ 
cies were reported only in the Tokachi region, Hokkaido; the possibility of its introduced origin from 
northeastern China is pointed out (see text). 
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Molecular phytogeny 

Leaves of the sample plants were dried and 
preserved in silica gel. After powdering the 
leaves in a mixer (Multi Bead Shocker MB455U, 
Yasui Kikai Corp, Osaka, Japan), total genomic 
DNA was extracted from 0.01 g of the powder us¬ 
ing a slightly modified version of the cetyltri- 


methyl ammonium bromide (CTAB) method de¬ 
scribed by Doyle & Doyle (1987). Polymerase 
chain reaction (PCR) was performed for the ITS 
and the trnK regions using the following primers 
[ITS: 1TS5 and ITS4 (White et al. 1990), trnK: 
trnKll and matK510R (Young et al. 1999)]. The 
PCR reaction mixtures contained 50-100 ng of 


Nil-Electronic Library Service 






The Japanese Society for Plant Systematics 


June 2014 


Fujii & AL.— Relict Distribution in Pedicularis Spicata 


template DNA, 6 pL of 10 x PCR buffer, 0.2 mM 
of each deoxyribonucleotide, 2.0 mM of MgCE, 
0.4 pM of each of the primer pairs, and 1.0 U of 
ExTaq DNA polymerase (Takara Bio Inc., Tokyo) 
in a total volume of 60 pL. The PCR program ran 
for 3 min at 94°C for initial denaturation, fol¬ 
lowed by 30 cycles of denaturation at 94°C for 1 
min, primer annealing at 53°C for 1 min and ex¬ 
tension at 72°C for 2 min. The reactions were then 
extended by 7 min at 72°C. PCR products for di¬ 
rect sequencing were excised from 1.5% agarose 
gels and purified using the GENECLEAN III Kit 
(Qbiogene, Inc., Carlsbad, CA). Sequencing reac¬ 
tions were performed using the CEQ dye termi¬ 
nator cycle sequencing Quick Start kit (Beckman 
Coulter, Inc., Fullerton, CA) according to the 
manufacturer’s recommendations. The following 
internal sequencing primers were used to read the 
whole regions: ITS: ITS3 and ITS4 (White et al. 
1990), trnK: matK-53F (Young et al. 1999), 
trnK600R (newly designed primer, 5’-AACTCG- 
GAGACAGGTAAACTT-3’) and matK840R 
(newly designed primer, 5’-TGTCATAACCTG- 
GATTTTCC-3’). The sequencing reaction prod¬ 
ucts were purified, concentrated by EtOH precip¬ 
itation, and then subjected to the CEQ 8000 Ge¬ 
netic Analysis System (Beckman Coulter, Inc.). 

The resulting nucleotide sequences were ed¬ 
ited using ChromasPro ver. 1.42 (www.technely- 
sium.com.au/ChromasPro.html) and aligned us¬ 
ing MEGA ver. 5.05 (Tamura et al. 2011). Phylo¬ 
genetic analysis was conducted using the neigh¬ 
bor-joining (NJ) method using MEGA ver. 5.05 


79 

and the maximum parsimony (MP) method using 
PAUP* ver. 4.0bl0 (Swofford 2003). In NJ analy¬ 
sis, genetic distances were calculated for both re¬ 
gions (ITS and trnK) using the Kimura two pa¬ 
rameter model (Kimura 1980) Insertions/dele¬ 
tions (indel) and missing data were treated as 
complete deletions. Clade support was estimated 
using bootstrap analysis (Felsenstein 1985) based 
on 10,000 replicates. In the MP analysis, the trees 
were solely constructed on the basis of substitu¬ 
tion data and using the branch and bound search 
option with MulTrees on. The addition sequence 
option was set at ‘furthest’ Clade support was es¬ 
timated using bootstrap analysis based on 1,000 
replicates in PAUP* ver. 4.0bl0. In the analysis, 
the heuristic search option was used, with the 
TBR branch swapping option and the ‘simple’ 
stepwise addition option. 

Morphological analyses 

Seven characteristics were measured on dried 
specimens to determine the taxonomic affiliation 
of the Hakusan populations: (a) angle of apex of 
the galea (upper lip of the corolla), (b) corolla 
length (excluding the calyx tube), (c) leaf blade 
length, (d) width between rachis and longest pin¬ 
na, (e) width between rachis and sinus of the pin¬ 
nae, (f) number of pinnae, and (g) plant length 
(Fig. 3). In addition, we measured the largest 
flowers, leaves, and plants. The measurements, 
excluding plant length, were taken using a porta¬ 
ble digital scope (Handy Scope 130s, Science- 
Eye, Saitama) with DinoCapture 2.0 ver. 1.3.5 



Fig. 3. Floral and leaf blade character¬ 
istics measured for analysis of 
variance, a: angle of the apex of 
galea, b: corolla length excluding 
calyx tube, c: leaf blade length, d: 
width between rachis and longest 
pinna, e: width between rachis 
and sinus of the pinna. In addi¬ 
tion to above-mentioned charac¬ 
teristics, we measured number of 
pinnae (f), plant length (g), and 
calculated the ratios of the char¬ 
acteristics c/d (h) and d/e (i). The 
flower and leaf in the figure are 
from plants of the Hakusan popu¬ 
lations (No. 1, Table 1). 
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Fig. 4. The neighbor-joining (NJ) tree using internal transcribed spacer (ITS) of nuclear ribosomal DNA among the species 
around Pedicularis spicata. Genetic distances were calculated on the basis of the Kimura two-parameter model. Numbers 
after the species or population names represent the individual numbers in Table 1. The maximum parsimony (MP) tree 
showed the same topology as the NJ tree. Numbers along the branches denote the bootstrap values (%); the former values 
are based on NJ analysis, while the latter values are based on MP analysis. 


software (AnMo Electronics Corp., Taipei). A 
measuring tape was used to determine the above 
ground length of the plants. We also calculated 
the ratio of characteristics using leaf characteris¬ 
tics (c-e) as follows: (h) ratio of leaf blade length 
to width between rachis and longest pinna, and (i) 
ratio of width between rachis and longest pinna to 
width between rachis and sinus of the pinnae. 
Characteristic h represents leaf shape, while 
characteristic i represents the degree of dissec¬ 
tion of the leaf pinnae. The software JMP ver. 
6.0.3 (SAS Institute Inc., Cary, NC, USA) was 
used for one-way analysis of variance (ANOVA) 
and Tukey-Kramer HSD tests were used to assess 
the significance of differences among species 
and/or populations. Furthermore, principal com¬ 
ponent analysis (PCA) was performed. Eigen val¬ 
ues greater than 0.5 were used as the criteria for 


extracting components. In our analysis, we used 
the data from characteristics a-g, while the ratio 
of characteristics (h and i) were excluded. 

We analyzed the morphological characteris¬ 
tics of the Hakusan populations using specimens 
and samples collected in the present study, except 
for characteristics of the capsules and seeds. 

Results 

DNA sequence variation within ingroup species 
Individual ITS sequences ranged in length 
from 571 to 575 bp. The aligned ITS data set was 
578 bp in length (including outgroups). The align¬ 
ment included the 5.8S ribosomal RNA coding 
region (164 bp). Pairwise distances between spe¬ 
cies/populations as determined using MEGA 
from the uncorrected /?-distance matrix ranged 
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Fig. 5. The neighbor-joining (NJ) tree using the trnK region of chloroplast DNA among the species around Pedicularis spicata. 
Genetic distances were calculated on the basis of the Kimura two-parameter model. Numbers after the species or popula¬ 
tions names represent individual numbers in Table 1. The maximum parsimony (MP) tree showed essentially the same 
topology as the NJ tree. Numbers along the branches denote the bootstrap values (%); the former values are based on NJ 
analysis, while the latter values are based on MP analysis. 


from 0.2% between accessions of Pedicularis 
verticillata (No. 6, 7 vs. No. 8) to 5.9% between 
accessions of P. verticillata (Nos. 6, 7, and 12) 
andP. refracta (Nos. 13 and 14) (Table 1). 

Sequences of the trnK region used in the pres¬ 
ent study ranged from 1026 to 1036 bp, and the 
aligned data set was 1044 bp in length. The align¬ 
ment included 377 bp of the matK coding region. 
The DNA sequences of the Hakusan populations 
were identical with Pedicularis spicata in this re¬ 
gion. Intraspecific variation in P. refracta was 
not detected. Values of the pairwise ^-distance 
ranged from 0.1% between accessions of P. verti¬ 
cillata (Nos. 7 and 8 vs. Nos. 9 and 11) and those 
ofP. verticillata andP. lineata (No. 11 vs. No. 16) 
to 1.9% between those ofP verticillata andP. re¬ 
fracta (Nos. 1-4 vs. No. 12). 

Phylogenetic analyses 

The NJ tree for ITS exhibited the following 
phylogenetic relationships (Fig. 4): First, the in¬ 
group accessions, including the Hakusan popula¬ 
tions, Pedicularis spicata, P. refracta, P. verticil¬ 
lata, P. lineata, and P. szetschuanica, formed a 


clade with a 99% bootstrap value. Second, two 
major clades were identified within the ingroup 
accessions. One clade comprised the Hakusan 
populations, P. spicata, and P. refracta with a 
97% bootstrap value. Within the clade, the Haku¬ 
san accessions and P. spicata were monophyletic 
with a bootstrap value of 94%. The other clade 
includedP. verticillata, P. lineata, andP. szetsch¬ 
uanica, supported by a bootstrap value of 44%. 
Within this clade, the accessions ofP. verticillata 
were monophyletic with 99% bootstrap support, 
and those of P. lineata and P. szetschuanica 
formed a clade with a bootstrap support of 100%. 
In the MP analysis for the ITS region, three MP 
trees, requiring 185 steps, were obtained. The 
consistency index, including uninformative char¬ 
acters (Cl), was 0.8486 and the retention index 
(RI) was 0.8627. The strict consensus tree of the 
three trees had the same topology as the NJ tree 
(Fig. 4). 

The NJ tree constructed using the trnK region 
of cpDNA is depicted in Fig. 5. The ingroup ac¬ 
cessions constituted a single clade with 96% 
bootstrap support. Bootstrap support for branch- 
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Table 2. Analysis of variance of morphological characteristics between Hakusan populations, Pedicularis spicata, P. verticil- 
lata, and P. refract a. 



Hakusan 

population 

P. spicata 

P. verticillata P. refracta 

a) angle of apex of the galea (°) 

52.9±12.5 C 1 
(N=20) 

50.8±10.6 C 
(N=27) 

83.9±11.3 A 
(N=95) 

66.9±10.7 B 
(N=23) 

b) corolla length (mm) 

8.9±0.8 C 
(N=20) 

10.1±1.8C 

(N=26) 

11.0±1.3 B 
(N=92) 

15.1±2.3 A 
(N=23) 

c) leaf blade length (mm) 

40.4±10.8 A 
(N=20) 

43.2±17.5 A 
(N=34) 

20.0±7.7 B 
(N=66) 

34.6±7.6 A 
(N=16) 

d) width between rachis and longest pinnae (mm) 

6.8±2.2 AB 
(N=20) 

5.4±1.8 BC 
(N=34) 

4.4±1.8C 

(N=66) 

0.8±0.4 C 
(N=66) 

e) width between rachis and sinus of the pinnae (mm) 

2.5±0.9 A 
(N=20) 

1.8±0.7 B 
(N=34) 

0.8±0.4 C 
(N=66) 

1.7±0.5 B 
(N=16) 

f) number of pinnae 

18.0±3.1 B 
(N=20) 

23.6±6.5 A 
(N=34) 

14.2±4.4 C 
(N=60) 

14.3±2.0 BC 
(N=15) 

g) plant length (cm) 

22.3±5.4 B 
(N=19) 

36.2±15.9 A 
(N=44) 

12.2±5.4 C 
(N=126) 

17.8±6.4 B 
(N=23) 

h) ratio of the characteristics (c) to (d) (c/d) 

6.2±1.4B 

(N=20) 

8.0±2.3 A 
(N=34) 

4.7±1.3 C 
(N=66) 

4.7±1.0 C 
(N=16) 

i) ratio of the characteristics (d) to (e) (d/e) 

2.8±0.6 C 
(N=20) 

3.2±1.3 C 
(N=34) 

5.9±1.9 A 
(N=66) 

4.8±1.3 B 
(N=16) 


1 Mean±S.D. (N=sample size), Letters (A—C) indicate significant differences at p < 0.05 determined by Tukey-Kramer HSD 
test. 


es within the ingroup was generally weak (48— 
69%), with the exception of a clade including Pe¬ 
dicularis verticillata and P. lineata (100%). The 
Hakusan populations and P. spicata were sister to 
this clade with 48% bootstrap support. In the MP 
analysis for the trnK region, 27 MP trees, requir¬ 
ing 90 steps, were obtained. Cl was 0.9222 and 
R1 was 0.8814. The strict consensus tree of the 27 
trees had essentially the same topology as the NJ 
tree, although several weakly supported branches 
in the NJ tree were unresolved polytomies (Fig. 

5). 

The topologies of the tree using the ITS re¬ 
gion and the trnK region were significantly dif¬ 
ferent from each other, particularly in Pedicu¬ 
laris lineata (Figs. 4 and 5). The incongruence 
length difference (ILD; Farris et al. 1994) test re¬ 
vealed that the ITS and the trnK region data were 
significantly incongruent (p = 0.001). 

Morphological analyses 

Using herbarium specimens and our field col¬ 
lections, we compared nine morphological char¬ 
acteristics among the Hakusan populations and 


the three species of Pedicularis (Table 2). The an¬ 
gle of the apex of the galea (a) did not differ sig¬ 
nificantly between the Hakusan and P. spicata 
samples (mean, 52.9°and 50.8°, respectively). 
Those values were smaller than for P. verticillata 
and P. refracta. The corolla (b) of the Hakusan 
samples was the shortest (mean, 8.9 mm) among 
the four samples, while the corolla of P. refracta 
was the longest (mean, 15.1 mm). The length of 
the leaf blade (c) of the Hakusan populations, P. 
spicata, and P. refracta was almost the same 
(mean, 34.6-43.2 mm), while the leaf blades of P. 
verticillata were the shortest (mean, 20.0 mm). 
The Hakusan specimens were relatively wider 
between the rachis and longest pinna (d) (mean, 
6.8 mm), but not significantly different from P. 
spicata. The width of the lamina between the ra¬ 
chis and sinus of the pinnae (e) was greatest 
(mean, 2.5 mm) in samples from Hakusan. The 
number of pinnae per leaf (f) was highest in P. 
spicata (mean, 23.6) and least in P. verticillata 
(mean, 14.2). In addition, P. spicata were the larg¬ 
est (mean, 36.2 cm), while P. verticillata was the 
smallest (mean, 12.2 cm). The ratio ofc to d, rep- 
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resenting leaf shape (h), the Hakusan populations 
and P. spicata showed relatively large values 
(mean, 6.2 and 8.0, respectively), while for i, the 
ratio of d to e, which represents the degree of dis¬ 
section of the leaf pinnae, the values for the 
Hakusan populations and P. spicata were rela¬ 
tively small (mean, 2.8 and 3.2, respectively). 

To determine the relationship of the Hakusan 
populations with Pedicularis spicata, P. verticil- 
lata and P. refracta, we performed PCA of the 
seven characteristics (a-g) mentioned above 
gathered from 104 specimens. A scatter diagram 
using principal components 1 and 2 is shown in 
Fig. 6. Component 1 explained 51.0% of the total 
variation. Leaf blade length (+0.49238 in factor 
loading) and the width between the rachis and si¬ 
nus of the pinna (+0.44382) contributed to this 
component (Table 3). Component 2 explained 
19.1% of the total variation, with the width be¬ 
tween the rachis and the longest pinna (+0.58291) 
and corolla length (+0.63099) contributing to this 
component. The cumulative contribution of the 
two components was 70.1%. Most specimens of 
P. verticillata, P. refracta, and P. spicata clus¬ 
tered within the range for each species, although 
specimens of P. verticillata and P. refracta were 
continuously distributed (Fig. 6). Specimens of 
the Hakusan populations were continuously dis¬ 
tributed with those of P. spicata. Furthermore, a 
gap was observed between the Hakusan popula¬ 
tions-/*. spicata cluster and the P. verticillata-P. 
refracta cluster. The samples from P. spicata 
from Hokkaido were continuous with the conti¬ 
nental populations of P. spicata. 

Description of the Hakusan populations 

Herbs, annual, erect, 10-30 cm tall. Caudex 
with simple slender root and one to several stems. 
Stems simple or sparsely branched, with four 
lines of hairs. Leaves often in whorls of four. Pet¬ 
iole of cauline leaves 0.5-1 cm long, narrowly 
winged, pubescent. Leaf blade oblong-lanceolate 
to narrowly lanceolate, 3-5 cm long, 0.5-1 cm 
wide, pinnately parted or pinnately lobed, apex 
acute; pinnae acute to cuspidate, with callous 
teeth, sparsely pilose on both surfaces. Inflores¬ 
cence terminal, spicate, 3-4 cm long at flowering 
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time, densely many flowered; proximal bracts 
leaf-like, distal bracts petiolate, lanceolate or 
rhomboid, acute, 5-10 mm long, 3-5 mm wide, 
callously cuspidate-dentate, pilose. Flowers July 
to August. Calyx suburceolate, ca. 4 mm long, 2 
mm wide, pilose, obscurely five-toothed. Corolla 
pink or purple, ca. 13 mm long; corolla tube ca. 9 
mm long; galea, ca. 4 mm long, 1.5 mm wide, 
apex rounded to truncate; lower lip ca. 5 mm 
long, 8 mm wide, trilobed, lobes orbiculate, ca. 4 
mm long, margin glabrous. Stamens 4; filaments 
4, 2 pubescent and 2 glabrous. 

Discussion 

Attributes of the Hakusan populations o/Pedicu- 
laris 

In the molecular phylogenetic analysis using 
ITS of nuclear ribosomal DNA (nrDNA), the 
Hakusan populations formed a strongly support¬ 
ed clade with Pedicularis spicata (Fig. 4). Al¬ 
though the Hakusan populations have been iden¬ 
tified as P. verticillata (Hakusan Nature Conser¬ 
vation Center 1995), all samples from Hakusan 
were positioned in another clade. P. refracta, P. 
lineata and/ 5 , szetschuanica were also positioned 
in a different clade. For the trnK region, the se¬ 
quences of the Hakusan populations and P. spi¬ 
cata were identical (Fig. 5). The results indicate 
that the Hakusan populations are nested within 
the P. spicata clade and should be considered to 
be members of that taxon. 

As mentioned above, the Hakusan popula¬ 
tions were identified as Pedicularis verticillata 
(Hakusan Nature Conservation Center 1995), but 
there are significant morphological differences 
between the Hakusan populations and P. verticil¬ 
lata in all nine characteristics examined using 
ANOVA (Table 2). In contrast, no significant dif¬ 
ferences were detected between the Hakusan 
populations and P. spicata with respect to the fol¬ 
lowing characteristics: (a) angle of the apex of the 
galea, (b) corolla length, (c) leaf blade length, and 
(d and i) width between rachis and longest pinna. 
According to Yamazaki (1993), P. spicata is dis¬ 
tinguished from P. verticillata and P. refracta as 
follows: in the corolla, apex of galea rounded or 
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Fig. 6. Scatter diagram showing 
principal components 1 and 2 
for 104 specimens of the 
Hakusan populations (□), Pe- 
dicularis spicata (■), P. verti- 
cillata (★), and P. refracta (A) 
(see Table 3). The two samples 
of P. spicata with arrow indi¬ 
cate that these localities are at 
Hokkaido, Japan. 


Compornent 1 (51.0%) 


truncate in Pedicularis spicata, but truncate in P. 
verticillata and P. refracta. The leaves are pin- 
nately lobed or parted in P. spicata, while pin- 
nately parted or dissected in P. verticillata and P. 
refracta. In our study, no evident differences in 
the apex of the galea were observed between the 
Hakusan populations and P. spicata. Further¬ 
more, no significant differences between them 
were observed in the degree of dissection of the 
leaf pinnae [i: ratio of the characteristics (d) to (e)] 
(Table 2). In addition, the results of PCA showed 
overlap in the morphological clusters represent¬ 
ing the Hakusan populations and P. spicata (Fig. 
6), whereas a much clearer gap was evident be¬ 
tween the Hakusan populations-/- 1 spicata clus¬ 
ter and the clusters representing P. verticillata 
and P. refracta. These findings show the Hakusan 
populations to be morphologically closer to P. 
spicata than to P. verticillata or P. refracta. 

Comparison of the description of the Hakusan 
populations and descriptions of Pedicularis spi¬ 
cata by Tsoong (1963), Yamazaki (1993) and 
Yang et al. (1998) did not reveal distinct morpho¬ 
logical differences. Based on the molecular and 


morphological evidence, we conclude that the 
Hakusan populations should be treated as P. spi¬ 
cata. 

Morphological differentiation o/Pedicularis spi¬ 
cata between Plakusan Mtns. and continental 
eastern Asia 

In external morphology, we observed signifi¬ 
cant differences in several characteristics be¬ 
tween the Hakusan plants of Pedicularis spicata 
and P. spicata from the Asian mainland (e, f, g, h; 
Table 2). The leaves of the Hakusan plants are 
smaller than in P. spicata from the mainland. The 
mainland plants have narrower leaves, but the 
Hakusan have fewer pinnae, although the length 
of the leaves was nearly the same in both. These 
results suggest a possibility of recognizing the 
Hakusan populations as an infraspecific taxon of 
P. spicata. Two subspecies of P. spicata were de¬ 
scribed by Tsoong (1963): subsp. bracteata and 
stenocarpa. The former is distinguished from 
subsp. spicata by its strongly enlarged bracts (10 
mm long x 8 mm wide), while the latter is distin¬ 
guished by the elongated capsules (8-10 mm 
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Table 3. Factor loading, eigenvalue, contribution, and accu¬ 
mulated contribution of principal components. 


Character 

Components 

PCI 

PC2 

a) angle of apex of the galea 

-0.382 

0.177 

b) corolla length 

-0.05150 

0.631 

c) leaf blade length 

0.492 

0.095 

d) width between rachis and longest 

0.339 

0.583 

pinnae 

e) width between rachis and sinus of 

0.444 

0.194 

the pinnae 
f) number of pinnae 

0.374 

-0.343 

g) plant length 

0.397 

-0.258 

Eigenvalue 

3.568 

1.334 

Contribution (%) 

50.976 

19.050 

Accumulated contribution (%) 

50.976 

70.976 


long) (Yang et al. 1998). Both the subspecies oc¬ 
cur in western Hebei, China (Tsoong 1963, Yang 
et al. 1998). The bracts are not enlarged in the 
specimens from the Hakusan populations (5-10 
mm long, 3-5 mm wide). However, we did not 
measure the attributes of the capsule. For the cap¬ 
sule, field research must be performed when the 
plants are in fruit. Thus, we were unable to de¬ 
finitively establish the intraspecific status of the 
Hakusan populations of P. spicata in the present 
study, and this needs to be addressed in a future 
study. 

Perspectives on disjunct distribution of Pedicu¬ 
laris spicata 

Pedicularis spicata occurs mainly in north¬ 
eastern Asia (Russia, Mongolia, China, and 
North Korea) (Fig. 2). In Japan, P. spicata has 
been reported only from the Tokachi district of 
Hokkaido. Yokoyama (1951) and Takita (2001) 
postulated that it had been introduced from north¬ 
eastern China, and it is treated as naturalized in 
several reports (Sugimoto 1978, Igarashi 2001). 
The evidence for its introduced status is based on 
only a few collections obtained from around a 
warhorse supplementary place and on roadsides 
in the vicinity of Ashoro and Nukabira, Tokachi 
district, Hokkaido [August 16, 1936, H. Yokoya¬ 
ma 3222 (TI, SAPS); September 17, 1983, K. 
Takita 5877; and August 21, 1994, K. Takita 5993 
(SAPT)]. We included these specimens in our 
morphological analyses, and most characteristics 


were within the range of variation of continental 
specimens (Fig. 6), but the evidence is inconclu¬ 
sive to determine its introduced status on Hok¬ 
kaido. 

This is the first report of the natural occur¬ 
rence of Pedicularis spicata in central Honshu, 
Japan, revealing a disjunct distribution pattern 
between the Asian mainland and Honshu Island. 
The studies of the biogeographic importance of 
the alpine region of central Honshu by Fujii & 
Senni (2006) and Ikeda et al. (2008) inferred that 
the alpine region of central Honshu acted as a re- 
fugium for alpine plants during the interglacial 
periods of the Pleistocene. Our study suggests 
that this may be the case for P. spicata. Based on 
the assumption that P. spicata originated on the 
Asian mainland, it is likely that the Hakusan pop¬ 
ulations represent surviving remnants of plants 
that reached Japan during cooler glacial periods. 

It is unclear why Pedicularis spicata occurs 
solely on the Hakusan population, while other 
high mountains in central Honshu also appears 
suitable. While many alpine plants occur on 
Hakusan Mtns., none are endemic (Shimizu 
1992, the Hakusan Nature Conservation Center 
1995), with the exception of Cirsium occidentali- 
nipponense Kadota on Mt. Sannomine in the 
southern part of the Hakusan Mtns. (Kadota 
1997) and some taxa at the rank of form (Satomi 
& Hashimoto 1970). The main alpine areas of 
Hakusan Mtns. are composed of volcanic rocks 
from eruptions that have taken place from ap¬ 
proximately 150,000 years ago to the present. 
The Hakusan populations of P. spicata, however, 
occur on relatively old volcanic rocks or on old 
sedimentary rocks of the Jurassic-Cretaceous 
Tetori Group (Geological Survey of Japan 2012). 
The Hakusan populations are found on gravel ex¬ 
posures, a habitat that other plant species are un¬ 
able to occupy easily. These factors may reflect 
the relict occurrence of the Hakusan populations 
of P. spicata, although they are insufficient for 
explaining why the species is absent elsewhere in 
Japan. 
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